Ground-State Entanglement in a Coupled-Cavity Model 
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Bipartite entanglement entropies are calculated for the ground state of the two-excitation sub- 
space in a two-site coupled cavity model. Each region in the phase diagram (atomic insulator, 
polaritonic insulator, photonic superfluid, and polaritonic superfluid) is found to be characterized 
by unique entanglement properties. In particular, the polaritonic superfluid region exhibits multi- 
partite entanglement among the two atoms and two cavity fields. This system provides a toy model 
in which a number of intriguing aspects of entanglement can be studied, such as the relationship of 
entanglement to phase transitions, entanglement of particles with different dimensionality, and the 
connection between experimentally accessible local observables and entanglement entropies. 
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Quantum phase transitions have been studied for a 
number of years in condensed matter physics [l[ . Unlike 
classical phase transitions, quantum phase transitions 
occur at zero temperature and are driven by quantum 
rather than thermal fluctuations. As a parameter is var- 
ied through the quantum critical point, the ground state 
of the system undergoes a transition characterized by a 
change in the correlations among particles. For exam- 
ple, the Bose-Hubbard model exhibits a transition from 
a product of localized states of definite particle number 
(Mott insulator) to a delocalized state with large fluctu- 
ations in the particle number per site (superfluid) Q. 

Examining the nature of these correlations led to the 
idea that entanglement could be a useful tool for under- 
standing the complex collective phenomena involved in 
quantum phase transitions d, 0, 0, H, 0, • This insight 
has been employed quite successfully in spin-^ systems, 
although other types of systems are more difficult to treat 
within the current theory of entanglement 0, Q . 

Recently it has been shown that an array of optical 
cavities, where each cavity contains one or more atoms, 
can un derg o a quantum phase transition 
3, Ell Ea E3, El El- Models of this tyP e can be ef- 



fectively described b y t he Bose-Hubbard model in some 
limits 1(| 11, 12, ID, 15 1 . However, the particles involved 
are not bosons but polaritons, quasiparticles composed of 
an atom plus photons. This difference adds a new dimen- 
sion to the problem 0, El, El El • 

Already there has been some work on the emergence of 
entanglement in coupled-cavity systems. Most of it has 
focused on the dynamic production of various entangled 
states for applications in quantum information process- 
ing. Proposals for creating entanglement between two 
polaritons 21 1, two photons 22, 23|, two atoms (23[, or 
a photon and an atom [23[ have appeared. There have 
also been proposals for creating highly ent ang led cluster 
states of polaritons 24, [5H|. Huo et al. [l5j compared 
the atom-atom entanglement and the visibility of photon 
interference fringes with the fluctuation of the total ex- 



citation number as means of characterizing the quantum 
phase transition. They found that the phase transition 
affects the entanglement between two atoms. 

In this paper I examine the ground-state entanglement 
in a simple two-site coupled-cavity model and how it re- 
lates to the various states in the phase diagram. Each 
site consists of a single-mode cavity containing a two- 
level atom; the two sites are identical. The cavities are 
coupled so that allows photons can hop between them. 
Although it is not strictly correct to refer to "phases" in 
such a small system, it has been shown that many results 
from finite models can be extrapolated to the thermody- 
namic limit [||. Crucially, in a system of this size it is 
possible to calculate the entanglement entropy of all the 
possible bipartitions of the ground state. Such detailed 
examination rapidly becomes more difficult as the size of 
the system increases. Furthermore, models with only a 
few sites could easily be implemented in superconducting 
qubit or quantum dot experiments 26|, 27, 



The Hamiltonian for the model is given by (h = 1) 



H = lu c a]aj + w„ \ej) {ej\ 

J'=l,2 

+ 9{a) \gj) (ej\ + a, \ej) + A(a\a 2 + a\a x ), 

(1) 

where lo c and uj a are the cavity and atom frequencies, g is 
the atom-cavity coupling strength, and A is the hopping 
strength. The operator dj (a -) is the lowering (raising) 
operator for the field in cavity j. The states \gj) and 
\ej) represent the ground and excited states of the atom 
in cavity j. Hence the operator \gj) (ej\ (\ej) (gj\) is the 
atomic lowering (raising) operator. 

A useful property of Eq. ([TJ) is that the total number 
of excitations N — a\di + d\a2 + |ei) (ei| + |ea) (e 2 | is 
conserved. For simplicity, the analysis is restricted to 
the case of two excitations. 

In the limit of zero hopping (A = 0), the system re- 
duces to a pair of noninteracting atom-cavity systems, 
each described by the Jaynes-Cummings model 2t| 3fl |. 
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The usual rotating- wave approximation has been applied, 
so that the interaction between the atom and the field 
contains only "energy-conserving" terms. This model is 
exactly solvable; the eigenstates of the individual cavities 
are given by 
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|0i) = \9i) \0i) 



(2) 



k > = sin y l e *) \ ( n - - cos y \9i) k) , (3) 

6 9 

\ n t) = cos y K n ~ + sin Y ' ^ 



(n = 1,2,3,...) is 
= 2gy / n/A, where 



where i = 1,2 denotes the site, |n) 
a photon number state, and tan 9 n 
A = io a — oj c is the detuning. 

The excited states are known as dressed states or po- 
laritons, and their nature depends on the atom-cavity 
detuning A. Consider, for example, the state At 
exact resonance (A = 0), it is a maximally entangled 
state of an atom and a single photon. As A — > +oo, 
| IT") — > |sf{) |1»): that is, the excitation becomes purely 
photonic. Likewise, as A — * — oo, the excitation becomes 
purely atomic. 

Previous work has shown that the coupled-cavity 
model can undergo a phase transition similar to that of 
the Bose-Hubbard model 



ii, m 



Near resonance 

the atom-cavity interaction shifts the frequency of the 

32] that 



cavity, causing the "photon blockade" effect [31 



prevents an additional photon from entering the cavity. 
The result is an effective on-site repulsion, as in the Bose- 
Hubbard model. For small values of A, this effective po- 
tential enforces the insulatorlike character of the ground 
state, characterized by a well-defined number of excita- 
tions at each site. As A increases, the hopping is able to 
overcome the repulsion. The ground state becomes delo- 
calized and the local excitation number is no longer fixed, 
similar to a superfluid state. Since only photons are able 
to hop between cavities, the limiting state as A becomes 
large may be characterized as a photonic superfluid. 

Figure [T] shows a phase diagram of the coupled cavity 
system. Four different regions have been identified, dis- 
tinguished by the phase (insulator or superfluid) and the 
nature of the particles involved [l(|. The insulator state 
may be polaritonic (for small A) or atomic (for large neg- 
ative A). As A becomes large, the ground state becomes 
a photonic superfluid. Perhaps the most interesting re- 
gion, however, appears in the vicinity of the transition 
around A = — A. The state in this region has both su- 
perfluid and polaritonic characteristics. 

The pha se diagram shown in Fig. Q] was obtained, as 
in Ref. [16j , by considering the variances of the total and 
atomic excitation numbers at a single site. For the Bose- 
Hubbard model, it has been shown that the entangle- 
ment properties of the ground state mirror the behav- 
ior of the number variance in signaling the phase transi- 
tion [tJ ■ Similarly, it seems reasonable to expect that the 
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FIG. 1: Phase diagram of the two-cavity system. The phase 
boundaries are defined as one-half the maximum value of the 
appropriate order parameter. For the superfluid state the 
variance of the total excitation number at site 1 is used as the 
order parameter, while for the polaritonic state the variance 
of the atomic excitation number at site 1 is used Hal. 



phases in the coupled-cavity model may be distinguished 
by their entanglement properties. 

There are two sources of entanglement in the cou- 
pled cavity system. The first is the polaritonic entangle- 
ment due to the atom-cavity interaction. From Eqs. ([3])- 
(TJJ it may be seen that the degree of entanglement of 
the polaritonic states is controlled by the interaction 
strength g and the detuning A. The polaritonic entan- 
glement is local in character in that it occurs between 
the atom and the cavity field at a single site only. Thus 
it is not included in entanglement measures that treat 
the interacting atom and field at one site as a single 
"particle". The second source of entanglement is the 
phase transition, controlled by the hopping parameter 
A. Phase-transition-induced entanglement may appear 
at any value of A, although the value of A at which the 
transition occurs depends on A. This type of entangle- 
ment is generally nonlocal in character. 

In this paper I focus on the bipartite entanglement in 
the ground state of the system, which is of course a pure 
state. Accordingly, the von Neumann entropy is a good 
entanglement measure. If a pure state of a composite sys- 
tem AB is given by the density matrix pab, the entropy 
S{pa) = — Tr(/5^ log 2 pa) of the reduced density matrix 
Pa = Ttb(pab) measures the entanglement between part 
A and part B of the system. 

Assuming that the two sites are identical, there are five 
possible bipartitions, as illustrated in Fig. [2] Each bipar- 
tition corresponds to an entropy measure. The single 
site entropy S(paici) is found by tracing over the atom 
and cavity of site 2 and applying the entropy formula 
to the resulting reduced density matrix. This entropy 
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— single site 

single atom 
— single cavity 

two atoms 

cross- site 



S(p,) 



FIG. 2: Illustration of the five different bipartitions in the 
two-site system. 



is essentially insensitive to the polaritonic entanglement. 
It therefore provides a measure of the nonlocal entan- 
glement. S(pai) and S(pci) give the entanglement of a 
single atom and a single cavity, respectively. We must 
also look at the entropy of both atoms S(paiA2) (equal 
to the entropy of both cavities) and the cross-site entropy 

S(pAlC2)- 

When comparing the various entropies, it is important 
to note that the von Neumann entropy has a maximum 
value of log 2 d in a <i-dimensional space. For a single 
two- level atom, d = 2 and S ma x{pAi) = 1- Although the 
dimension of a single cavity field is infinite, the restriction 
to two excitations means that each field forms an effective 
three-level system, so S max (pci) = log 2 3 ~ 1.585. The 
atomic subspace has dimension 4, giving S ma x{pAiA2) — 
2. The single-site and cross-site bipartitions each contain 
one atom and one cavity field for a nominal dimension of 
6. However, the restriction to two excitations limits the 
number of atom-cavity states to 5, which corresponds to 

SmaxiPAlCl) = S max (pAlC2) = log 2 5 1=3 2.322. 

Figure [3] shows the five bipartite entropies in the limit 
of small hopping. At large negative detuning all of the 
entropies are nearly zero, as expected from the charac- 
terization of this region as an atomic insulator. Near 
A = 0, four of the five entropies reach a local maxi- 
mum, reflecting the polaritonic entanglement created by 
the resonance between the atom and the cavity field at 
each site. However, the single site entropy remains close 
to zero and thus there is no intersite entanglement. This 
corresponds to a polaritonic insulator state, which takes 
the form of a tensor product across sites where each site 
contains an entangled state. Finally, as A increases, the 
single atom and atomic subspace entropies drop to zero 
while the other three entropies asymptotically approach 
1.5. In this photonic superfluid region, the atoms are 
both in the ground state. The photons in the two sites are 
entangled, as evidenced by the nonzero single site, single 
cavity, and cross site entropies (which become equivalent 
as A — > oo). 

The term "polaritonic superfluid" suggests a delocal- 
ized state composed of polaritons, which are themselves 
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FIG. 3: Bipartite entropies as a function of A in the small- 
hopping limit, A = O.Olg. 
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FIG. 4: Bipartite entropies as a function of A in the large- 
hopping limit, A = lOg. 



entangled states of an atom and a cavity field. This sug- 
gests that the polaritonic superfluid state might exhibit 
true multipartite entanglement. Although a multipartite 
entanglement measure for four particles of mixed dimen- 
sionality does not yet exist to our knowledge, an indicator 
of multipartite entanglement is that all possible biparti- 
tions of the system are entangled. A plot of the entropies 
in the large-hopping limit is shown in Fig. [4j As in the 
small- hopping limit, the state at large negative detuning 
is an atomic insulator, while at large positive detuning 
the state is a photonic superfluid. However, all five bi- 
partite entropies reach a maximum near the transition 
point —A w A, indicating that there is indeed multipar- 
tite entanglement in the polaritonic superfluid state. 

Note that although only two extreme cases typifying 
the small- and large-hopping limits are shown, the en- 
tanglement peaks follow the transition throughout the 
whole parameter space. At intermediate values of A, the 
entropy curves transform between the two limiting cases 
in a predictable fashion. 

It should also be pointed out that the multipartite en- 
tanglement in this system only occurs because the atom- 
cavity structure allows multiple subsystems at each lat- 
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tice site. Other authors have considered the relation be- 
tween phase transitions and multipartite entanglement 
among more than two lattice points 0, [H, 34 , 3{| , but 
that is a somewhat different problem. 

Entanglement properties, therefore, provide a useful 
means of differentiating the various regions in the phase 
diagram. The atomic insulator state is completely sepa- 
rable. In the polaritonic insulator state, polaritonic en- 
tanglement exists between the atom and field at each 
site, but there is no intersite entanglement. The photonic 
superfluid state exhibits bipartite entanglement between 
the sites. Genuine multipartite entanglement appears 
in the polaritonic superfluid state. This characteriza- 
tion suggests a connection between the entropies and the 
number variances that have served as order parameters 
in previous work, which is currently being investigated. 
However, it seems that no single bipartite entropy pro- 
vides a complete view of the system; they need to be 
taken together for a full description. 

Furthermore, the description in terms of entanglement 
reveals a qualitative difference between the transitions 
at large and small hopping values. In the small hopping 
limit, the change in particle type from atomic to polari- 
tonic to photonic is due to the resonance between the 
atoms and photons. The phase transition occurs at a 
larger value of A than the resonance, so the region of po- 
laritonic behavior does not overlap much with the super- 
fluid region. As a consequence, the entanglement remains 
primarily bipartite. In contrast, for large hopping values 
the phase transition occurs deep in the negative detun- 
ing regime. The change in particle type is driven by the 
transition itself, resulting in multipartite entanglement 
in the transition region. Whether this difference, and in- 
deed the polaritonic superfluid state itself [la |. persists 
in the thermodynamic limit remains an open question. 

The model studied here, although itself quite simple, 
opens up some significant avenues of inquiry. It provides 
a toy model for the problem of entanglement of parti- 
cles with different dimensionality The interplay between 
on-site and intersite entanglement is another interesting 
problem: how is entanglement created, transformed, and 
destroyed in the system? Of course, there are also ques- 
tions about practical applications, particularly whether 
this model can be used to produce useful entanglement 
for tasks such as computation. Furthermore, the connec- 
tion between the variances and entropies could poten- 
tially offer an experimentally accessible way to measure 
entanglement [361 ] . 
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